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Abstract. We propose a new experimental test bed that uses ions in the 
collective ground state of a static trap to study the analogue of quantum- 
field effects in cosmological spacetimes, including the Gibbons-Hawking effect 
for a single detector in de Sitter spacetime, as well as the possibility of 
modeling inflationary structure formation and the entanglement signature of 
de Sitter spacetime. To date, proposals for using trapped ions in analogue 
gravity experiments have simulated the effect of gravity on the field modes by 
directly manipulating the ions' motion. In contrast, by associating laboratory 
time with conformal time in the simulated universe, we can encode the full 
effect of curvature in the modulation of the laser used to couple the ions' 
vibrational motion and electronic states. This model simplifies the experimental 
requirements for modeling the analogue of an expanding universe using trapped 
ions, and it enlarges the validity of the ion-trap analogy to a wide range of 
interesting cases. 
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1. Introduction 

Curved-spacetime quantum field theory (QFT) is an essential ingredient in our modern 
understanding of cosmology [1] and has also given rise to remarkable insights and puzzles 
concerning the nature of Planck-scale physics [2]. Unfortunately, most predictions of curved- 
spacetime QFT are extremely difficult to test directly, confined as we are to a low-curvature 
region of spacetime, while dimensional analysis conspires to make the most straightforward 
effects extremely tiny. This situation has given rise to a number of proposals for testing analogue 
curved-spacetime QFT effects, where spacetime is replaced by a system of coupled atoms, 
and field modes are replaced by the collective normal modes of oscillation of the atoms — 
effectively a 'phonon field'. The advantage is that while the observable effects are conceptually 
the same (analogous), the parameters of the experiment may be adjusted to result in a much 
stronger and more easily observed effect. The most common analogue curved-spacetime QFT 
proposals involve Bose-Einstein condensates (BECs), liquid helium or ion traps as the analogue 
spacetime [3, 4]. 

Here, we employ a general but often overlooked concept to extend the range of possible 
analogue spacetime QFT experiments in ion traps, and we discuss several analogue cosmology 
experiments that appear to be greatly simplified by the new picture, proposing one of them 
(simulation of the Gibbons-Hawking effect [5]) in detail. The main idea, which has been 
proposed previously in the context of BEC simulations [6, 7], is to encode the effects of 
spacetime curvature into a modulation of the analogue detectors, rather than in the motion 
of the atoms that comprise the analogue spacetime. We will show how this amounts to 
making a shift in the analogy away from laboratory time -o- detector proper time and towards 
laboratory time conformal time. In general, such a shift in analogy can be made between the 
experimental lab time and any new time coordinate in the simulated spacetime, but in the case 
of Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmologies in the conformal vacuum, 
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we will see that this particular choice (laboratory time made analogous to conformal time) 
encodes the entire effect of the spacetime curvature and therefore enables us to propose analogue 
curved-spacetime QFT experiments in which the analogue spacetime is fixed (resembling flat 
Minkowski space), thus making our proposals far more accessible experimentally. 

Building on efforts to implement elementary quantum gates, ion trapping now offers some 
of the most exquisite demonstrations of quantum control in modern physics. The ability to 
couple the internal electronic state of one or more trapped ions to the collective vibrational 
degrees of freedom enables the ability to produce a large class of pure quantum states [8]-[10], 
including Fock states, squeezed states and cat states, as well as quite complex information 
processing tasks [11, 12]. We use the control of the trapping potential afforded by ion traps, 
together with the ability to reach quantum limited motion, to propose the construction of 
new ion-trap analogues of the single-detector response in FLRW cosmology (including the 
Gibbons-Hawking effect [5]), along with several proposals for ways to apply the model to 
situations involving multiple detectors, including inflationary de Sitter cosmology [13, 14] 
and the entanglement signature of de Sitter spacetime [15]. We discuss the most important 
experimental parameters governing the simulation and the regime in which the analogue 
experiment is likely to conform closely to the simulated cosmology. Highly efficient quantum 
measurements, based on fluorescent shelving [8], provide a practical means to test our 
predictions. The main advantage of trapped-ion analogue gravity experiments over other 
proposals, such as BECs [6, 7], lies in the exquisite quantum control over state preparation, 
interaction and readout. In particular, the fact that each ion is equipped with its own highly 
efficient and tunable motional detector (through laser coupling [10]), which counts individual 
phonons, means that ion traps have a particular advantage in efficient quantum control of the 
detector interaction and readout [3]. We compare our proposal to previous work [16, 17] that 
suggests using ion traps with a time-dependent trapping potential to simulate cosmological 
particle creation. 

The paper is organized as follows. Section 2 describes in more detail the conceptual shift 
from the standard 'field picture' of analogue QFT to the 'detector picture', which we use 
extensively, and reviews the relevant properties of FLRW cosmology that make it possible. 
In section 3, we review the fundamental features of ion-trap simulations of cosmology, which are 
essential for describing our proposals. Section 4 is devoted to the description of our experimental 
proposals. Our concluding remarks are contained in section 5. Unless otherwise specified, we 
use natural units, in which Ti = c = 1. 

2. Conformal invariance: field picture and detector picture 

To build an experimental analogy to quantum field effects in spacetime, we will make use of 
two basic ingredients: (i) a coupled, discrete system of oscillators, where deviation in a given 
oscillator from its equilibrium position q„(t) plays the role of the field 0(x, t), and (ii) one or 
more detectors, which can be coupled to the analogue field in a way that models a given type of 
particle detector in spacetime. In the case of the ion traps that we discuss in the next section, for 
example, the electronic state of a specific atom is coupled to the deviation from its equilibrium 
position. 

Before we discuss the analogue experiment, however, we recall a few basic properties of 
standard Unruh-De Witt detectors in FLRW spacetime [18]. Recall that the usual coupling of 
the field to detector is expressed via the following interaction Hamiltonian (which, importantly, 

New Journal of Physics 12 (2010) 095019 (http://www.njp.org/) 



4 



IOP Institute of Physics <J> Deutsche physikalische gesellschaft 



assumes that the cosmic time variable t equals the proper time r for the detector), 
Hmt(t) = f(t)m(t)<p(x,t), 



(1) 



where all operators are in the interaction picture, m(t) is the detector monopole moment and 
/(f) is a 'window function' modeling a time-dependent coupling of the field to the detector. 
In the case of a single detector with energy gap £ at a fixed location x, this leads to a detector 
response function given by 



This is the probability that the detector, starting from its ground state \g), will become excited 
into its higher-energy state \e) after interacting with the field on application of the Hamiltonian 
H int (t) over a very long time — i.e. at least long enough to cover the (nontrivial) support of 
fit). This formula assumes that the interaction is weak enough for time-dependent perturbation 
theory to be used. Here, D + (t, t') is the two-point function given by 



for a given state of the field (recall that the detector is fixed at spatial position x), and the bar is 
used with foresight to indicate that the time coordinate being used is t (cf equation (7) below). 

An analogue experiment is one that reproduces the essential features of the response 
function in curved spacetime for an analogous detector coupling and the state of the oscillators 
(typically with the ground state of the collective normal modes of the oscillators being analogous 
to a vacuum state of the spacetime field). What we would like to emphasize here is that this can 
be achieved in multiple ways: by manipulating the collection of oscillators themselves or by 
manipulating the detector. When the spacetime involved is conformal to Minkowski spacetime, 
as in FLRW cosmology, the entire effect of curvature can be encoded completely within either 
the oscillators or the detector, leading to two distinct 'pictures' for analogue experiments. 

2.1. The field picture 

In the field picture, the relevant features of the analogue spacetime are encoded in the oscillators, 
and thus directly in the modes of the analogue field and the resulting two-point function. 
Thus, in the response function, the analogy is contained in the two-point function directly, 
i.e. (0(x, t)4>(x, t')) <h>- {(p„(t)(p n (t')) (where 0„(f) is an appropriately dimensionless version of 
q„(t)). The analogy is completed by physically moving the oscillators — an expanding array of 
ions simulating an expanding spacetime, for example [16, 17]. The important feature of this 
picture is that t labels both the proper time of the analogue detector in the lab and the proper 
time of the detector being simulated. 

The field picture is thus immediate conceptually, and for this reason it has been the 'default' 
picture for proposed analogue cosmological QFT experiments in the literature. We emphasize, 
however, that it is not the only picture for cosmological models and that adequate control of the 
analogue spacetime (the oscillators) often presents a set of significant experimental challenges. 

2.2. The detector picture 

When a given spacetime is conformal to Minkowski spacetime, and when the fields considered 
are massless and conformally coupled, we can switch to the detector picture, in which the array 




(2) 



D + (f,r') = (0(x,O^(x, t')) 



(3) 
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of oscillators is held fixed (in analogy to Minkowski spacetime), but the detector is modified in a 
time-dependent fashion to create the same analogue effect as a curved spacetime. The important 
conceptual feature of this picture is that the laboratory proper time describing the evolution of 
the analogue detector is not analogous to the proper time of the simulated spacetime detector. 
Instead, the experimental proper time is analogous to the conformal time of the simulated 
spacetime and detector. This is accomplished by introducing a time-dependent energy gap in the 
analogue detector, which simulates an ordinary, constant-energy gap of a detector in a curved 
cosmological spacetime. The payoff for the conceptual shift and the modulation of the detector 
is a greatly simplified treatment of the analogue spacetime — the array of oscillators is held fixed. 

To see how this works, we recall the transformation properties of fields in spacetime 
(see [18] for more details). For simplicity, we work with a scalar field. Recall that a spacetime 
(with metric g llv ) conformal to Minkowski spacetime (wih metric h /lv ) can be coordinatized 
such that g^vix, £) = £2(x, t;) 2 h llv . (In this section x = (x, £), with £ playing the role of a time 
coordinate.) In these coordinates, the two-point function in the conformal vacuum satisfies 

D + [x(0, f ; f '1 = ^(x) (2 -" )/2 L»;[x(0, f ; *(?'). tW) {2 - n) i\ (4) 

where Z)j^[x(f), £ ; x(f '), f] is the standard Minkowski two-point function for a scalar field 
along the trajectory x(£), n is the number of spacetime dimensions, and we omit the bar over 
D + because these are functions of £ instead of t (cf equation (2) above and equation (7) below). 
The variable £ is thus interpreted as the conformal time of an FLRW universe on the left-hand 
side of the above equation and the Minkowski time on the right-hand side. 

In FLRW cosmology with no spatial curvature, the metric in comoving coordinates takes 
the form 

n-1 

di 2 = df 2 -a(0 2 ^(dx') 2 , (5) 

i=i 

where t is cosmic time and also the proper time for a comoving observer (i.e., an observer on 
an inertial trajectory of constant x). Transforming to conformal coordinates by substitution of 
the conformal time £ , which satisfies 

dt 

d£ = — , (6) 
a(t) 



results in a metric of the following form, 



ds 2 = £2(£) 2 



d? 2 -£(d*') 2 



(7) 



where 



£2(0 =a[f(?)]. (8) 



By using equation (4), the two-point function on the trajectory of a comoving detector at a fixed 
location x can be expressed as 

^FLRWC* 0 = £>FLRW(^ fO = « (O^^fi (O^'^C?, O. (9) 

where £> + (£, £') is short for £> + [x(f ), ?; x(f ), £'] for both the FLRW and Minkowski Green's 
functions as functions of £, and we explicitly show the connection to the Green's function 
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D + (t,t') used in equation (2). The factor of f2(f ) (2 -" )/2 £2(f ') (2 ~" )/2 may be absorbed into 
the definition of the detector window functions, and thus an analogue Minkowski two-point 
function is equally good for describing FLRW cosmology, provided we can re-interpret the ex- 
perimenter's 'clock time' as the conformal time f . For this interpretation to work, however, the 
detector must be designed to evolve with respect to the simulated proper time t as e~ lE ' in order 
to model a comoving detector, rather than evolving in the simulated conformal time f as e~ l£? . 

The appropriate modification of the detector requires that we simulate time evolution in 
t with a t -independent Hamiltonian H 0 by means of a ^-dependent Hamiltonian i/o(f )■ This 
time-dependent Hamiltonian will be the physical Hamiltonian acting on an individual oscillator 
in the lab [6, 7, 19]. To find H 0 (t;), we simply write the desired Schrodinger equation in t and 
transform to f . Thus, a time-independent H 0 in the equation 

d 

i— * = (10) 
df 

can be re-expressed as 

i-jU = # 0 (f )*, (11) 
df 

where H 0 (t;) = (df/df )H 0 . Simulating a detector with fixed energy gap E in t, given by H 0 , 
requires a time-dependent 'laboratory Hamiltonian' H 0 (t;), which corresponds to a simple time- 
dependent scaling of the energy gap in laboratory time. For the case of an FLRW cosmology 
simulation, this is given by 

H Q (O = V(t;)H Q = a[t(OlH 0 . (12) 

Observe how these ingredients come together to describe the response of a comoving 
detector at a fixed location x in an FLRW universe, 

At / d^cO/CrV^'-^o^LRwC*. t') 

-OO J — CO 

= / df/ df'^^ 7 /[ f (0]/[f'(f')]e- l£[f(f) -''«' )] 



'?(/=-oo) Ji;(t=-ao) 

x £2(f) (2 -" )/2 £2(0 (2 "" ,/2 £>M(^ O 



/■?(f=oo) fi(t=oo) 

= df / df 'F(OF (f ')e-™-''«' )] D^ , O, 

Jt(t=-oo) Jt(t=-oo) 



• f(r=oo) /•f(r=cx>) 

(13) 

'f(/=-oo) Jf(/=-cx>) 

where 

F(f ) = ^C2(0 (2 -" V 7[*(f)] = a[f (? )] (4 "" )/2 /[/(f )] (14) 
df 

is the window function in f (conformal time, also Minkowski time) necessary to simulate the 
window function f(t) in cosmic time. Note the non-standard exponential factor in the last line, 
coming from the f -dependent energy gap. This connection between the field picture and the 
detector picture is illustrated schematically in figure 1 . 

To summarize, the behavior of a detector with fixed energy gap E, determined by H 0 , 
on a comoving trajectory in an FLRW universe in the conformal vacuum may be simulated 
by a detector with a time-dependent energy gap, determined by H 0 (t;), moving inertially in the 
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(a) Field picture (b) Detector picture 

Figure 1. Schematic representation of the field picture and detector picture for 
a detector in an expanding universe exhibiting conformal invariance. In the 
field picture (a), expansion of the universe (shown left to right) causes the field 
modes to redshift, while the detector resonant frequency remains unchanged. 
The field modes are defined with respect to cosmic time t, using the metric from 
equation (5), and the detector Hamiltonian H 0 is independent of t, satisfying 
equation (10). Conformal invariance allows an equivalent description in the 
detector picture (b), in which expansion of the universe (shown left to right) 
causes the detector resonant frequency to blueshift, while the field modes remain 
unchanged. The field modes in this picture are defined with respect to conformal 
time £, using the metric from equation (7), and the detector Hamiltonian H 0 (t;), 
given by equation (12), is now f -dependent, satisfying equation (11). With 
additional modulation of the detector window function (equation (14)), both 
pictures result in exactly the same prediction for the physical response of the 
detector, as shown in equation (13). 



Minkowski vacuum. In addition, the detector's window function must also be modified: /(f) m>- 
F(t;). The relevant equations for these substitutions are equations (10) and (14), respectively. 
This correspondence applies equally well to multi-detector effects, such as entanglement, and it 
may be exploited in the construction of analogue gravity experiments. 

3. Cosmology in an ion trap 

With the power of the detector picture, a vast number of gravitational models become 
experimentally accessible simply by modulating the detector parameters appropriately. Our 
analogue gravitational model of choice is the linear ion trap [10]. Previous authors have 
considered ion traps for simulating de Sitter spacetime [16] (including the Gibbons-Hawking 
effect [5]) and cosmological particle creation due to a finite-duration period of expansion [17]. 
We will compare our proposal with these in due course. We also note that trapped ions 
rotating around a ring trap have been proposed for detecting the acoustic analogue of 
Hawking radiation [20]. In what follows, we describe the ion-trap analogues of the field 
(in the Minkowski vacuum), an ordinary De Witt monopole detector [21] and a modulated 
detector. 

We will go into some detail in describing the analogy between trapped ions and fields 
in spacetime in order to make certain nuances clear, which can have an important effect 
when generalizing the model to backgrounds other than Minkowski spacetime. The reader is 
encouraged to note, specifically, that the actual analogue of the field 0 is a dimensionless discrete 
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field 4>„ (equation (19)); that the usual way of writing the ion-laser interaction (equation (22)) 
conflates interaction-picture time dependence (equation (20)) with additional time dependence 
from the laser, which has important consequences when trying to apply equation (10); that the 
Lamb-Dicke parameter rj (equation (23)) becomes time dependent when the trap frequency 
changes 4 and when the laser frequency changes (equation (29)); and that changing the laser 
frequency may, in some cases, also affect the Rabi frequency S2 0 . again through equation (29). 
We restore explicit use of ti in this section. 



3.1. Discrete phonon field 

In a linear ion trap, a string of up to about ten ions is confined to one dimension using 
a combination of radio-frequency and static electric fields and capped at the ends by static 
electrodes [8]. Linearizing the total potential created by the harmonic potential due to the trap 
electrodes and the mutual Coulomb repulsion between the ions, we can model the collective 
motion of the ions in a linear trap as a collection of coupled harmonic oscillators. With a 
coupling matrix A obtained from the linearized combination of these potentials, the Hamiltonian 
is given by [10] 

1 _,rj- . Mv 2 „ 

H= p T p + ^A^, (15) 

2M l 1 2 H 1 

where q = {q\ q N ) T is a column vector of operators corresponding to the displacement of 

each ion from its equilibrium position, p = (pi, . . . , p N ) T is a column vector of corresponding 
momenta, M is the mass of each ion and v is the effective harmonic trap frequency provided 
by the trap electrodes (typically, v is approximately a few MHz [8]). The coupling matrix 
A is symmetric and positive-definite. Thus, we can diagonalize it as A = B T AB, where A = 
diag(/Lii, . . . , fi N ) is a diagonal matrix of positive eigenvalues, arranged in ascending order. 
The orthogonal matrix B defines the transformation to normal-mode coordinates Q = Bq and 
momenta P = Bp. Equation (15) can be rewritten in these coordinates as 

N P 2 Mv 2 N 

P =\ P =i 

where 

v P ■= v^fjJ^ (17) 
is the oscillation frequency of normal mode p and (since the normal modes oscillate 
independently) we have diagonalized this Hamiltonian using the standard prescription for a set 
of independent harmonic oscillators, ignoring the zero-point energy. The normal-mode raising 
and lowering operators satisfy [a p , a p ,] = 8 PP <, as is appropriate for independent bosonic modes. 
The local oscillations q„, about the ions' equilibrium positions are given, in the interaction 
picture, by the following equation [10], 

N i: 

- iv ^a p + e w ^a\), (18) 

where m 6 {1, . . . , N] labels the ion and is an entry of B, which defines the spatial mode 
functions of the normal modes 5 . We can define a unitless version of this local displacement 

4 While this does not affect our proposal, it does affect the papers [16, 17] and any other field-picture proposals 
for ion-trap simulations of gravitational effects. 

Prefactor and phase conventions for q„, vary in the literature. 
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operator as well, 

JL b ( P> 

p=\ 11 p 

with 5 m (f) = *J~h/2Mv(p m (t;). This discrete field will serve as the analogue of the Minkowski 
scalar field <p(x, f ). Note that laboratory clocks tick off conformal time, so the time dependence 
of these modes is as it would be for modes in Minkowski spacetime, and the time coordinate is 
appropriately labeled by £ . Because the coupling matrix A only approximates a simple 'balls on 
a string' coupling (characteristic of a bosonic field), the field will have a non-standard dispersion 
relation, given by the frequencies v p from equation (17), as a function of the normal-mode 
index p. This deviation from ideality means that the analogy with a massless scalar field is an 
imperfect one. 



3.2. Laser-induced coupling of vibrational and electronic states 

A feature of this model is that the internal electronic state of each ion can be turned into a local 
detector for the displacement of that ion. To achieve this for a given ion, its internal state must 
become correlated with its linear displacement from equilibrium. This correlation is provided 
by an external laser, which is used to drive an electronic transition between two meta-stable 
electronic levels \g) and \e), separated in energy by Tlcoa [8, 10, 22]. 

The interaction between an external classical laser field and the mth ion is described in the 
dipole and rotating- wave approximation 5 by the interaction-picture Hamiltonian [8, 10], 



Hf'iO = -ihQo )e lMl * C0, »*' l »- fl »-f J - a™ {S )e 



(m) ( r\p~^ k COS 6 1m (f )-<»Lf ] 



(20) 



where fi 0 is the Rabi frequency for the laser-atom interaction (typically, fi 0 ~ 100 kHz [8]), 
k is the magnitude of the laser's wave vector k, which makes an angle 9 with the trap axis, 
q m (£) is the interaction-picture position operator for the mth ion, and o± (?) = e ±I<BAf o± 
are its interaction-picture electronic raising and lowering operators, where cr+" > = \e) mm {g\ and 
crl" 1 ' = \g) mm {e\. We may additionally absorb the laser's time-dependent phase e T1£0Lf into the 
interaction-picture operators £± m) (? ) := e^ 1 ^ cr^^) = e ±lA? cr± , where 

A = co A - co L (21) 

is the detuning of the laser below the atomic transition. This allows us to write the interaction 
Hamiltonian as 

H\ m) (0 = -ihQo [ ? {'")(f )e +ltcos ^» tt) - ^")(^) e -'* cos ^»(f)J , (22) 
which is the common form in the literature [8, 10]. 

6 There are actually two different rotating-wave approximations that will be mentioned in this paper. The first, 
discussed only here, is between the laser frequency and the electronic levels of the ion. This one is always valid 
since the laser frequency &>l ~ co a S> v. In [8, 10], this is discussed in more detail. The second rotating-wave 
approximation — which is what will be meant by all further instances of this term — is between the detuning 
in equation (21) and the normal-mode frequencies. The validity of this approximation depends on these two 
frequencies, the coupling strength and the time of the interaction. See [23, 24] and references cited therein for 
more details. 
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The size of the rms fluctuation in q,„ as compared to the wavelength of the laser is measured 
by the Lamb-Dicke parameter, 



hk 2 cos 2 6 Ax£l 

V-J— T7} -^(27rcos0), (23) 

2Mv h 



where y/h/2Mv is the ground-state rms fluctuation of the center-of-mass mode of the ions. This 
quantity is representative of the overall rms fluctuations of the with ion, denoted by Ajc^, since 
the frequencies v p for all higher normal modes of the ions remain within an order of magnitude 
of v for small numbers of ions (up to N ~ 10), realistic for current experiments [10]. Typical 
values of the Lamb-Dicke parameter are r\ ~ 0.01-0.1 [25]. When n^Jn + 1 <3C 1, the so-called 
'Lamb-Dicke limit' (n is the average phonon number), the ion is well localized with respect to 
the wavelength of the laser, and we can expand the exponentials in equation (22) to first order 
in r], which is equivalent to first order in kcos 8 q m (f ), giving 

H\ m) (0 ~ fi Q 0 ff w (O + fi £2 0 fc cos 6 q m (?) g ( x m) (?) , (24) 

carrier sideband 

where = e +lA ^i"" +e- lAf CT l" ,) and g^(0 = -ie +lAf CT j m) + i e" lA ^o-l m) . The first term 

corresponds to excitation of the transition directly by the laser, while the second couples the 
atomic transition to vibrational motion. 

With the laser sufficiently detuned away from resonance (| A | > v) and having a sufficiently 
narrow bandwidth (interaction time longer than v~ [ ), we can neglect the carrier term and keep 
only the sideband term, 

H W (f) mQY] y „ (e . iVr(ar +e i„ rf t)^ w +e -^al m) ). (25) 

This Hamiltonian can be used to couple the atomic transition to the vibrational motion, with the 
detuning A serving as the effective resonant frequency for the detector. When A = v p (called 
the first red sideband), this interaction models a detector that becomes excited when a phonon 
is absorbed — an ordinary phonon detector. When A = — v p (called the first blue sideband), it 
becomes excited when a phonon is emitted — a rather unorthodox detector indeed. After these 
interactions have completed, readout (projection) of the electronic state is accomplished using 
a bright laser that strongly couples the excited state to an auxiliary level, causing excited ions to 
fluoresce, while ions in the ground state remain dark [10]. 

As long as the coupling strength Q 0 r] is weak enough, no more than a single \g) — > \e) 
transition will be excited for a given ion in the time during which the laser coupling is active. 
This allows us to use perturbation theory to calculate the detector response function. Strictly 
speaking, an infinite interaction time is needed for the simplified detector description given 
above to be exact. For finite times, the non-resonant terms (such as ajcr + ) can be non-negligible, 
resulting in red-sideband detections even when the ions are in the collective ground state. But 
this is what we want: this interaction has the form of a De Witt monopole coupling [21], and 
keeping these non-resonant terms is important for modeling the response of the modulated 
detector used to simulate the effects of expansion, the details of which we now discuss. 

3.3. Time-dependent laser coupling 

The usual construction, described above, produces a De Witt monopole detector for the ions' 
vibrational quanta. For our purposes, however, we need to introduce time dependence into our 
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detector. What we need, specifically, is a detector whose Hamiltonian is scaled by a[t(t;)]. But 
we cannot simply apply this directly to the laser interaction Hamiltonian because the effective 
energy gap of the 'detector' is fi A, which is the difference between the atomic transition energy 
and the laser energy. Therefore what we need is to apply the following modulation, 

Aw.A(?):=a[*(f)]A, (26) 

which, if we wish to vary only the laser frequency, can be achieved by 

fflLW-a)L(0:=a)A-fl[*(f)]A. (27) 

Alternatively, one could choose to Stark shift the atomic frequency such that co A i-> &>a(0 := 
&>l + a A. This may be more experimentally viable for some applications, but we will 
assume from now on that all modulation is done within the laser. 

The only other change necessary to simulate an FLRW metric in the detector picture is a 
new window function. Window functions correspond to modulation of the Rabi frequency Q Q 
of the atom, which is proportional to the amplitude of the interaction laser [10]. Since a window 
function was needed anyway for initiating and ceasing the laser interaction, the requirement of a 
new window function presents no additional fundamental challenge. These modifications result 
in the following interaction Hamiltonian, 

flf'(0 S ' debmd > hOoriFU) Y %(e- iVrf flr + e 1 ^a, t )(e lA ' (f Vf J +e- lA ' (f) al m) ). (28) 

One other relevant thing to note is that the Lamb-Dicke parameter r\ depends on the laser 
wavenumber k, which is proportional to o>l(0- Thus, modulating the laser frequency will also 
modulate the coupling strength automatically through a time-dependent r](£), 

— = ^ = ^^ = 1-— (o[t(f)]-l), (29) 
r] k co L q)l 

but since a> L ~ a> A ~ 1 GHz, while A ~ v ~ 1 MHz, the fractional change in the coupling 
strength due to the frequency shift will be ~ 0.1% unless the scale factor gets to be the inverse 
of that same order. This is unlikely, however, because such a large modulation (beyond a 
few tens of MHz) is both experimentally challenging [26, 27] and unnecessary because the 
detector would be so far detuned from the normal mode frequencies that it would (at best) 
detect nothing and (at worst) begin to excite extraneous atomic transitions. Therefore, while 
the time dependence of the Lamb-Dicke parameter must be considered in actual experiments, 
it can easily be compensated by suitably modulating the laser power. We also note that a similar 
argument may apply to the Rabi frequency £2 0 , depending on whether the atomic transition 
being excited is dipole-allowed or dipole-forbidden. In the case of the latter, Q Q also varies with 
the atomic transition frequency [10], which is of the same order as the laser frequency, and thus 
the calculation in equation (29) applies to Q (£) as well. Regardless of the experimental methods 
needed to achieve it, the new interaction Hamiltonian remains as in equation (28). 

At this point, we should take some time to talk about the approximations that could 
break down because of this additional modulation. One of the biggest concerns is avoiding 
unwanted atomic transitions. If the bandwidth of the modulation, which is governed by the 
inverse of the detection time (width of the window function), is too large, this broadband 
response could overlap with unwanted transition energies. In addition, the range over which 
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the frequency varies could cause the same problem, even if the modulation is slow and over 
a long period of time. This issue strongly depends on the spectrum of the actual atoms used. 
Typically [8], however, v ~ 1 MHz, while <w A ~ 1 GHz, which should provide quite a large 
bandwidth over which the laser frequency may be swept. Finally, the rate of modulation 
affects the feasibility of the scheme as well. Modulations up to about 20 MHz can be achieved 
with an acousto-optic modulator, while faster rates (up to a few GHz) can be achieved using 
an electro-optic modulator [26, 27]. We now describe several possible applications of ourFLRW 
analogue system. 

4. Experimental applications 

4.1. Two distinct mechanisms by which expansion creates particles 

Birrell and Davies [18] go to great lengths to distinguish between two effects that might be 
called 'particle creation' due to expansion. The first is the non-zero detection rate of an ordinary 
detector in what should be considered an 'empty' universe when the detector is on during the 
expansion. If the expansion eventually stops, and the universe is asymptotically Minkowski, 
this effect is no longer seen. This does not necessarily mean that particles are not detected in 
that case but just that this particular mechanism for particle production only applies during 
expansion. 

There is another type of 'particle creation' that might be observed after a finite period of 
expansion has passed. Imagine a universe starting asymptotically in the Minkowski vacuum for 
very early times. Now imagine that expansion begins and ends smoothly, leaving the system 
in a non-equivalent vacuum but in a spacetime that is also asymptotically Minkowski for late 
times. Particle creation occurs in this case because of the Bogoliubov transformation mixing 
annihilation and creation operators of the original vacuum into operators of a single type in 
the new vacuum (a a fter = a<2 be fore + ^a before ), tnus ensuring that there will be real particles flying 
around after the expansion. 

As we have formulated it, the detector picture can model the first case (detection during 
expansion) but not the second (detection after expansion). The reason is the need for conformal 
invariance 7 . We have explicitly required that the vacuum we model in the detector picture be the 
conformal vacuum, which is (by definition) the ordinary Minkowski vacuum in conformal time. 
If we were to model a finite-time expansion in the detector picture, particles would be detected 
during the expansion, but none would be detected afterwards 8 . 

An ion-trap analogy has been proposed for each of these effects. Alsing etal [16] attempted 
to model the thermal response of a detector in de Sitter spacetime (exponential expansion) by 
modulating the trap potential and thereby modifying the phonons' normal modes. Since the 
proposal involved simulating the expansion by modifying the dynamics of the ions, the proposal 
uses the field picture. There were several problems with the proposal, including incorrectly 
predicted dynamics, calculational mistakes and ions that would fly out of the trap in a realistic 

7 Our work does not rule out a detector picture for cases that do not admit conformal invariance. Our present 
formulation of it, however, requires this invariance. 

8 To see this, just imagine a detector whose resonant frequency is first fixed at one value, then smoothly shifted 
to another value and then fixed at that value, but only after this process is detector activated. For long enough 
detection times (much longer than the inverse detector frequency), the rotating-wave approximation may be made, 
and no particles will be detected [23, 24], 
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experiment. These problems were identified and discussed in [17] (see endnote [5] therein), 
which itself was an alternative proposal for witnessing particle creation in an ion trap — again 
by chirping the trap frequency — but this time the proposal was for the second type of particle 
creation in which real particles are created by the simulated expansion and detected after it 
has ceased. This proposal does not suffer from the conceptual, calculational and experimental 
problems of [16], but it also does not model the non-zero response of a detector during 
expansion. 

With the power of the detector picture, however, we can model a detector that is on during 
expansion — the original goal of [16] — but with a completely static trap, instead embodying 
the entire effect of the expansion within the modulation of the laser used to couple the ions' 
vibrational motion to their electronic states. In fact, we are not limited to de Sitter spacetime 
in this picture; we can model any scale factor as long as the required modulation is feasible 
experimentally. 



4.2. Single-detector response 



The detector response function A m for the mth ion, to lowest non-trivial order in the Dyson 
series expansion, is given by 



/oo /»oc 
dt / dC'F(f) J F(f')e- lA[ '« ) - t «' )] (^ m «)0m(?')). 
■OO J — OO 



(30) 



This is the probability that an ion will get flipped to its excited electronic state by interaction with 
the field. The interaction is assumed to be weak enough to only excite one transition from \g) to 
\e) during the detection, thus allowing the use of time-dependent perturbation theory [28]. We 
have extended the limits to ±oo, with finite limits able to be enforced using the window function 
F(t;), as required. Also note that the last term (in angular brackets) is just the ion analogue of 
Z)Jj(f , £')• Since the field is the ions' collective ground state, we can plug in equation (19) and 
evaluate 

N \b^f 

(31) 

-1 y/^P 



.(f)*»(?'))=E^L e " i, '"' c ' 



giving 



n— 1 V^P J-OO 



iAt(f) p-ivpt 



(32) 



Note that we can change integration variables back to t to obtain 

dt f(t) e 



P =i 



f 

J — c 



At e -iv p f (0 



(33) 



This formula assumes n = 2, which is the case for a (1 + 1) -dimensional spacetime — this being 
the most natural spacetime to simulate in a linear ion trap, which has one effective spatial 
dimension. We note, however, that one-dimensional effects in spacetimes equipped with more 
than one spatial dimension can be simulated as well, at the expense of additional powers of a (t) 
being included (see equation (14)). The second form, equation (33), may be interpreted as the 
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field-picture calculation of the same quantity as in the detector-picture formula, equation (32), 
since in this integral the detector time dependence appears ordinary, being just e~ lA ', while 
the normal modes now experience time-dependent modulation, as e -1 "''^". Also note that in 
the trivial case of no simulated expansion, f = t, and this quantity gives zero in the limit 
where f(t)^ 1 (unless A < 0, indicating a blue-sideband detector), as a result of the rotating- 
wave approximation being applied to a detector in the vacuum. We can, however, still achieve 
a non-zero response in this case for short enough detection times [23, 24]. 



4.3. De Sitter spacetime and the Gibbons— Hawking effect 

One of the most celebrated features of de Sitter spacetime is the thermal response of a local 
detector to the conformal vacuum with a temperature proportional to the expansion rate, 

7gh = ^, (34) 

in natural units. This is the Gibbons-Hawking effect [5]. It means that a single inertial detector 
in the conformal vacuum of de Sitter spacetime will respond as if it were bathed in thermal 
radiation at a temperature Ton in Minkowski spacetime. An important physical difference 
between these two cases is that in the case of Gibbons-Hawking radiation, every inertial detector 
believes itself to be at rest with respect to the thermal bath it responds to — i.e. no velocity- 
dependent Doppler shift is observed [5] — while a thermal state in Minkowski spacetime has a 
preferred rest frame, and detectors in relative motion with respect to this frame (even if still 
inertial) will perceive a Doppler-shifted spectrum [18]. 

It is this thermal response that [16] wanted to simulate by chirping the trap potential and 
modifying the phonons' normal modes. Since the proposal involved simulating the expansion by 
modifying the dynamics of the ions, it uses the field picture. Later work corrected this proposal 
for a single ion [29], which again involved exponentially modulating the trap frequency. The 
results were different from those calculated (incorrectly) in [16] but still showed an effectively 
thermal response with a perceived thermal state with temperature 7gh 9 - Here we will instead 
propose to model the Gibbons-Hawking effect in the detector picture. 

De Sitter spacetime is an FLRW universe with an exponentially increasing scale 
factor in t, 

a(r) = e" a[t(0] = -—, (35) 

with transformations to and from conformal time given, respectively, by 

f(0 = t(0 = (36) 

K K 

which means that small intervals in each variable are given by 

df=e""dr dt = - — , (37) 

where we have chosen t = 0 to correspond to £ = — 1 /k . (In addition, t — >• oo corresponds to 
£ -> 0~ and t —> —oo corresponds to £ —> —oo.) As in all simulations we are proposing, the 
clock on the wall in the laboratory ticks off seconds of conformal time £, while the simulated 

9 This is measured by taking the ratio of the detector response functions at the first red and the first blue sideband, 
giving e- 2 " A /" ) which is the desired ratio. 
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proper time x of the detector is equal to the simulated cosmic time t and is given by equation (36) 
in terms of £ . Ideally we would like our detector to have an infinite detection time in t to get a 
perfectly thermal spectrum. However, this is not physical, so we will simulate it with a detection 
time that is sufficiently long. Following the procedure proposed in [29] (but modulating the 
detector instead of the trap frequency), we will model the physical detector as having a flat 
window function over a finite detection time from 0 to T, 

f(t)~®(t)-®(t-T), (38) 

where @(m) = 1 for positive u and 0 for negative u. The equality is not strict because we want 
smooth turn-on and turn-off for the detector. The details of these functions are not important 
because the approximations we use allow us to model this detector as being on for a very long 
time 10 , 

v » k and v e"" T k. (39) 

A look at equation (33) reveals that in the field picture, the time dependence of the field modes 
in an FLRW universe satisfies q-' v p^'\ These limits, respectively, correspond to an expansion 
that begins in a regime where the rotating-wave approximation is valid (at t = 0) but extends 
into a regime where the field modes have effectively stopped oscillating (at t = T), meaning that 
particles will be detected [23, 24]. The first approximation lets the lower limit of the window 
function be extended from t = 0 to t — >• — oo because the detector will see nothing during f < 0 
anyway. The second allows the upper limit to be extended from t = T to t — >• oo because the 
detector, which is assumed to be tuned to a frequency A ~ v, is far off -resonance from the 
normal-mode frequencies v p , which are all within an order of magnitude of v [10]. 

In the detector picture, this corresponds to sweeping the trap frequency from A (the actual 
detector frequency that we want to simulate) at f = — 1/k to Ae* r at £ = —er KT /K, with 
intermediate values determined by 

A(?) = -A. (40) 

This is a polynomial modulation in f, while the proposals in the papers [29] and [16] 
require exponential modulation of the trap frequency. The window function in conformal time 
(assuming n = 2) is 

F(S) = ~- (41) 
Plugging these into the integral of equation (32) gives 

/o /.O 
df F(f)e- lA ' (f) e- iy " c = / df {-kO { ^' k) ~ 1 z~ iM ■ (42) 
-CO J — DC 

Substituting u = —k^, du = — /cdf, P = A/k and a = v p /ic, along with appropriate 
convergence factors (a — >• a +ie, /? — >• f$ — ie, with c — > 0 + ), gives 

Kj 0 ic(—iayp Ka ,p 

10 We use the approximations from [29] rather than those of [16] because those of the latter (which correspond to 
rapid expansion of the trap initially) do not actually allow the limits on the window function to be extended as the 
authors claimed. 
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Plugging in for the substituted variables and plugging the final result back into equation (32), 
along with the identity |r(i/6)| 2 = 7r/[/Jsinh(7r^)], gives 

A_(A)= — —-nr, L_^_L (44) 

kA e 2 * A /«- 1 ^fJT„ 

p=l v 1 

which is a Planck spectrum at the Gibbons-Hawking temperature. This formula may be found 
in [16], but it does not correspond to the procedure that the authors proposed. It is also different 
from the result obtained in [29] for the case of a single ion. This is expected, however, since 
modulation of the trap frequency creates single-mode squeezing that has no analogue in de 
Sitter inflation. 

There are several experimental hurdles to be overcome. Firstly, we want k to be low so 
that we can use the approximations in equation (39), but we do not want it to be so low that 
detection events are prohibitively rare (<1%). Secondly, the laser modulation rates must be 
within feasible limits. Thirdly, the electronic structure of the ions must be such that extraneous 
transitions are not excited while modulating the laser. We can address the first two problems 
briefly here while noting that any real application will require a more detailed analysis involving 
the atomic spectra. 

Combining the approximations from equation (39) gives er" T <JC 1, which means that 
q-kt ^ o.05 or less. This means that k^(T) ~ 0.05/cf (0) and so f (T) — £(0) ~ k~\ meaning 
that the detector has to be on for roughly the inverse of the expansion rate. This makes the 
bandwidth of the modulated laser spectrum ~k <3C v, so the normal mode frequencies are still 
well resolved, and we do not have to worry about broadband excitation of extraneous atomic 
transitions. The detector energy gap A(f ) (and thus the laser frequency <w L ) must be modulated 
over a range (e KT — 1)A ~ 20 MHz for A ~ v ~ 1 MHz [8]. The laser power similarly needs to 
increase by a factor of ~20 over the course of the simulation 1 1 . 

While all of the above is achievable with current technology [26, 27], the detection 
probabilities are prohibitively low due to the exponential scaling with A/k, which we have 
required to be large in equation (39) (since A ~ v). Nevertheless, having established the 
theoretical connection to Gibbons-Hawking radiation for always-on detectors in our model, 
higher detection probabilities can be achieved using window functions of shorter duration 
and/or by using larger values of k. The resulting response function will no longer be a Planck 
spectrum but can be calculated numerically. This situation still models the analogous de Sitter 
detector — just one with a finite-duration window function. 

There are two things we want, then, in an experimentally accessible simulation. The first is 
to make the absolute excitation probability (for the red sideband) high enough to be detectable; 
the second is to have some sort of signature of thermality. To boost the detection rate, we will 
have to violate the approximations from equations (39). This means that the limits of the integral 
in equation (42) will be finite, 



/ d£ f(Oe- lA ' (f) e-"'" f = - / duu 

Jt, K Ju, 



r(i/3, — iauf) — T(i/3, — \au{) 



tea 1 ? 



where = —Kt;^ = e and T(z, b) = f, dx x z 1 e 1 is the upper incomplete gamma 
function, which is a built-in special function in many computer algebra packages. We define the 

11 Simulating a higher-dimensional spacetime will change this requirement (see equation (14)). 
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time interval in terms of the simulated initial and final cosmic times — t Y and ff, respectively — in 
order to emphasize that this deviation from a perfect Planck spectrum is not an experimental 
artifact but results rather from simulating a detector that is on only for a finite duration 
(specifically, from t\ to tf) in de Sitter spacetime. Defining the regularized upper incomplete 
gamma function Q(z,b) := T(z, b)/ T{z), we can see that equation (44) gets modified to 

A m (A; tl ,t f )= y °" R(A; tl ,t f ), (46) 

kA e 2 * A / K - 1 



where 



^ \ K K J \ 



a: 



2 

, (47) 



which can be evaluated numerically. Extending the limits to ±oo gives A m {A) = 
A,„(A; —oo, oo), in agreement with equation (44). 

In the case of effectively infinite-duration detectors, we can take the ratio of the red- and 
blue-sideband detection probabilities and find the signature of thermality, A m (A) / A m (— A) = 
e -2irA/K p Qr fj n j te detection times that are shorter than required for the detector to 

thermalize with the perceived Gibbons-Hawking bath — which is the regime we are required 
to work in — this relation no longer holds. However, a new ratio can be calculated from 
equation (46), 

A,„(A;f„f f ) = e _ 2]lA/K R(A; t u t t ) 
A m (-A; tl ,t t ) R(-A;t u tf)' 

which depends on when the detector is switched on and off but once again agrees with the fully 
thermalized result in the limit of long detection time. 



4.4. Multiple detectors 

The ion-trap analogy using the detector picture carries over to the case of many simultaneous 
detectors without modification. This is because all terms in the Dyson series are integrals 
involving the interaction Hamiltonian, and such integrals can always be recast from cosmic 
time to conformal time, resulting in a new interpretation as the response of Minkowski- 
spacetime detectors with time-dependent frequencies, as long as the window function is also 
appropriately modified. This leads us to suggest the possibility of studying the two-point 
correlations of scalar-field modes in de Sitter spacetime by measuring the correlation function 
for simultaneous electronic excitations in multiple ions [28] under the detector modulation 
prescription proposed here. Such correlations are believed, in inflationary cosmology, to be 
responsible for structure formation in the early universe [13, 14]. Two-component BECs 
have served in recent work [31] as analogue systems for studying inflationary cosmology. 
Because of the restriction to conformally invariant settings in the current proposal, our goals 
are comparatively modest. By measuring correlation functions in the detection probability for 
excitations in multiple ions [28] under a detector modulation corresponding to various scale 
factors, we can analyze the spectrum of fluctuations that should result from the corresponding 
cosmological predictions [13, 14]. 

Independently of the cosmological implications, Retzker et al [32] propose a method of 
swapping ground-state entanglement of a chain of ions to the electronic states of the ions at the 
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ends of the trap using sonically local interactions (i.e. detection events that happen faster than 
the speed of sound can propagate across the trap). This is the ion-trap analogue experiment of 
the proposal in [33] to use the Minkowski vacuum of a quantum field to generate entanglement 
between distant local detectors through entirely local interactions with the field. Using the 
detector picture, we can extend the proposal of [32] — as well as any other proposal that seeks to 
model only the Minkowski vacuum — to an equivalent proposal that instead models any FLRW 
universe in the conformal vacuum. This is particularly exciting because it brings into the realm 
of possibility performing the analogue of the experiment proposed in [15] to distinguish heating 
of a Minkowski field from de Sitter expansion by comparing whether entanglement is generated 
between two local detectors. The detector interactions used in [15] are impractical both in the 
actual vacuum of outer space and in an ion trap because the detection probabilities would be 
far too low. However, in an attempt to implement this analogue experiment, the non-standard 
detection techniques of [32], which help boost the entanglement significantly over that achieved 
by the ordinary monopole detectors, may be applied to the following two cases: (a) modulated 
energy levels and window functions, corresponding to de Sitter expansion, and (b) a heated field 
with no detector modulation. We reserve the details of this proposal for a future paper, but we 
include it in the discussion to illustrate the wide ranging applicability of the detector picture in 
ion-trap simulations of quantum fields. 

5. Conclusion 

We have demonstrated an alternative 'picture' for analogue QFT in curved spacetime, which 
encodes the effect of curvature in the detector rather than in the motion of the analogue 
quantum field, and we have used this technique to propose new and more accessible analogue 
experiments for curved-spacetime QFT effects. These include the Gibbons-Hawking effect for 
a single detector in de Sitter spacetime, plus the possibility of modeling structure formation in 
inflationary cosmology and the entanglement signature of de Sitter spacetime. 

A question that may arise in the mind of the reader is the degree to which the ion 
trap/curved spacetime analogy continues to hold when time has been re-defined in such a 
way as to make the trap resemble Minkowski spacetime. As one example, ordinary Minkowski 
spacetime has no cosmological horizon — one of the most conceptually important ingredients of, 
for example, de Sitter spacetime. Careful consideration of points like this shows that the analogy 
still holds and that the horizon is actually still present. The horizon in the detector picture is 
related to the fact that evolution from — oo to +00 in de Sitter cosmic time corresponds to an 
evolution from —00 to 0 in the analogue Minkowski time variable £ . The cosmological horizon 
is thus re-interpreted as the backward light cone extending from the detector location at £ = 0 
in Minkowski spacetime, which equally well shields the detector from some of the field modes 
that are necessary to define the vacuum state. 

In a more general sense, however, all analogue experiments are, at their root, quantum 
simulations. In a simulated system, there is no reason to expect that simulated time must 
evolve in the same way as ordinary time. Anyone who has done numerical simulations of time- 
dependent systems can attest to this. Going between the two pictures can be thought of as 
running exactly the same simulation program on two different machines — while the simulation 
may be run at different rates (depending on the machine), the end result is precisely the same 
output, and the degree to which the output is analogous to something in the 'real world' is 
identical. 

New Journal of Physics 12 (2010) 095019 (http://www.njp.org/) 



19 

Acknowledgments 



IOP Institute of Physics <J> Deutsche physikalische gesellschaft 



We thank Greg Ver Steeg for stimulating discussions and ideas, Dave Kielpinski for consultation 
on the experimental feasibility of our scheme and Tim Ralph for a forum in which to present and 
discuss our ideas. We also thank the Defence Science and Technology Organisation (DSTO) 
for support, as well as the University of Queensland and the Australian Research Council 
Centre of Excellence for Quantum Computer Technology, whose hospitality towards NCM was 
instrumental in bringing this work to fruition. Research at the Perimeter Institute is supported 
by the Government of Canada through Industry Canada and by the Province of Ontario through 
the Ministry of Research and Innovation. 

References 

[1] Weinberg S 2008 Cosmology (New York: Oxford University Press) 

[2] Hawking S 1 975 Particle creation by black holes Commun. Math. Phys. 43 1 99 

[3] Schutzhold R 2009 Recreating fundamental effects in the laboratory? Adv. Sci. Lett. 2 121 

[4] Barcelo C, Liberati S and Visser M 2005 Analogue gravity Living Rev. Relat. 8 12 

[5] Gibbons G W and Hawking S W 1977 Cosmological event horizons, thermodynamics, and particle creation 
Phys. Rev. D 15 2738 

[6] Fedichev P O and Fischer U R 2003 Gibbons-Hawking effect in the sonic de Sitter space-time of an expanding 

Bose-Einstein-condensed gas Phys. Rev. Lett. 91 240407 
[7] Fedichev P O and Fischer U R 2004 Observer dependence for the phonon content of the sound field living on 

the effective curved space-time background of a Bose-Einstein condensate Phys. Rev. D 69 06402 1 
[8] Leibfried D et al 2003 Quantum dynamics of single trapped ions Rev. Mod. Phys. 75 281 
[9] Monroe C et al 1995 Resolved-sideband Raman cooling of a bound atom to the 3D zero-point energy Phys. 

Rev. Lett. 75 4011 

[10] James D F V 1998 Quantum dynamics of cold trapped ions with application to quantum computation Appl. 
Phys. B 66 181 

[11] Leibfried D et al 2003 Experimental demonstration of a robust, high-fidelity geometric two ion-qubit phase 
gate Nature 422 412 

[12] Schmidt-Kaler F et al 2003 Realization of the Cirac-Zoller controlled-NOT quantum gate Nature 422 408 
[13] Bardeen J M, Steinhardt P J and Turner M S 1983 Spontaneous creation of almost scale-free density 

perturbations in an inflationary universe Phys. Rev. D 28 679 
[14] Bassett B A, Tsujikawa S and Wands D 2006 Inflation dynamics and reheating Rev. Mod. Phys. 78 537 
[15] Ver Steeg G and Menicucci N C 2009 Entangling power of an expanding universe Phys. Rev. D 79 044027 
[16] Alsing P M, Dowling J P and Milburn G J 2005 Ion trap simulations of quantum fields in an expanding 

universe Phys. Rev. Lett. 94 220401 
[17] Schutzhold R et al 2007 Analogue of cosmological particle creation in an ion trap Phys. Rev. Lett. 99 201301 
[18] Birrell N D and Davies P C W 1982 Quantum Field Thoery in Curved Space (Cambridge: Cambridge 

University Press) 

[19] Olson S J and Ralph T C 2010 Entanglement between the future and past in the quantum vacuum 
arXiv: 1003.0720 [quant-ph] 

[20] Horstmann B et al 2010 Hawking radiation from an acoustic black hole on an ion ring Phys. Rev. Lett. 104 
250403 

[21] DeWitt B S 1979 Quantum gravity: the new synthesis General Relativity, An Einstein Centenary Survey 

(Cambridge: Cambridge University Press) pp 680-745 
[22] Wallentowitz S and Vogel W 1996 Motional quantum states of a trapped ion: measurement and its back action 

Phys. Rev. A 54 3322 

New Journal of Physics 12 (2010) 095019 (http://www.njp.org/) 



20 IOP Institute of Physics <J> Deutsche physikalische gesellschaft 

[23] Drummond P D 1987 Unifying the p ■ A and d ■ E interactions in photodetector theory Phys. Rev. A 35 4253 
[24] Zaheer K and Zubairy M S 1988 Atom-field interaction without the rotating-wave approximation: a path- 
integral approach Phys. Rev. A 37 1628 
[25] Walls D F and Milburn G J 2008 Quantum Optics 2nd edn (Berlin: Springer) 
[26] Salech B E A and Teich M C 1991 Fundamentals of Photonics (New York: Wiley) 
[27] Yariv A 1989 Quantum Electronics (New York: Wiley) 

[28] Menicucci N C and Milburn G J 2008 Spatial correlation functions for the collective degrees of freedom of 

many trapped ions arXiv:0807.3500vl [quant-ph] 
[29] Menicucci N C and Milburn G J 2007 Single trapped ion as a time-dependent harmonic oscillator Phys. Rev. 

A 76 052105 

[30] Turchette Q A et al 2000 Heating of trapped ions from the quantum ground state Phys. Rev. A 61 063418 
[3 1] Fischer U R and Schiitzhold R 2004 Quantum simulation of cosmic inflation in two-component Bose-Einstein 

condensates Phys. Rev. A 70 063615 
[32] Retzker A, Cirac J I and Reznik B 2005 Detecting vacuum entanglement in a linear ion trap Phys. Rev. Lett. 

94 050504 

[33] Reznik B, Retzker A and Silman J 2005 Violating Bell's inequalities in vacuum Phys. Rev. A 71 042104 



New Journal of Physics 12 (2010) 095019 (http://www.njp.org/) 



